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1.  Introduction 


The  Gruneisen  function,  T,  is  a  key  thermodynamic  function  in  the  study  of  materials 
subjected  to  the  high  pressures  of  shock  loading.  Macroscopically,  it  relates  the  pressure 
changes  to  changes  in  internal  energy  at  constant  volume,  by  way  of  the  relation 

r  =  v(dp/dE)v  ,  (i) 


where  p  is  pressure,  E  is  internal  energy,  and  V  is  the  volume  of  the  material's  unit  cell. 
In  general,  T  may  be  a  function  of  the  lattice  volume  (i.e.,  spacing  of  the  atomic 
oscillators)  as  well  as  the  temperature  (i.e.,  amplitude  of  the  vibration).  However, 
Gruneisen  theory  assumes  that  the  vibrational  amplitude  is  small  relative  to  the  lattice 
spacing;  thus,  the  set  of  oscillators  comprising  the  lattice  possess  frequencies  that  are 
assumed  to  vary  with  volume  alone.  In  terms  of  this  characteristic  vibrational  frequency 
of  the  lattice,  u,  the  volume  dependence  of  the  Gruneisen  function  is  expressed, 
according  to  classical  theory  (1,  2),  as 


F  jv  =  du/dV 

L 0 


(2) 


During  the  decades  in  which  the  classical  approach  to  lattice  vibration  held  sway,  a 
number  of  analytical  Gruneisen  formulations  were  presented  in  the  literature,  with 
several  gaining  a  widely  cited  reputation  (2,  3,  4).  Because  of  the  manner  in  which  these 
particular  models  were  formulated,  T  was  given  in  terms  of  both  the  material's  volume 
and  cold  pressure,  pc  (i.e.,  the  0°  isothermal  pressure).  However,  since  the  cold  pressure 
is,  itself,  a  function  of  volume  alone,  these  formulations  for  T  indeed  satisfy  the 
Gruneisen  assumption  of  temperature  independence. 

While  these  important  models  retained  the  key  definitional  coupling  embodied  in 
equation  2,  many  later  developers  of  numerical  hydrocodes  for  impact  simulation  opted 
to  bypass  the  formalism  of  these  historical  models,  in  favor  of  direct  empirical 
characterization  of  T(H)  according  to  preferred  functional  forms.  Perhaps  the  motivation 
in  this  was  that  empirical  shock  data  provided  a  material's  shock  Hugoniot,  and  not  the 
cold-pressure  curve  demanded  of  the  classical  models.  Regardless  of  the  motivation, 
once  analytically  decoupled  from  the  mechanics  of  lattice  compression,  these  wholly 
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empirical  formulations  for  T  were  not  always  thermodynamically  robust.  Segletes  (5-7) 
has  shown  how  these  decoupled  Griineisen  formulations  are  susceptible  to  an  unusual 
host  of  thermodynamic  inconsistencies  that  can  lead  to  nonphysical  thermodynamic 
instabilities  in  shock/impact  calculations. 

There  is  value,  then,  in  reconsidering  the  classical  Griineisen  methods  as  an  important 
bridge  between  modern  ab  initio  methods  and  the  decoupled  empirical  approaches  that 
still  characterize  many  extensively  utilized  numerical  codes  in  the  shock/impact 
community. 


2.  Classical  Griineisen  Theories 


The  classical  models  by  Slater  (2),  Dugdale  and  MacDonald  (3),  as  well  as  the 
tree-volume  theory  cited  by  Vashchenko  and  Zubarev  (4),  have  been  widely  studied  and 
utilized  and  relate  a  materials'  Griineisen  behavior  to  the  compression  behavior  of  the 
lattice.  While  it  is  somewhat  difficult  to  compare  the  original  derivations  of  these  three 
models  side  by  side  to  ascertain  the  particular  points  of  departure  that  bring  about  their 
distinctions,  some  have  attempted  to  do  exactly  this. 


While  the  forms  originally  presented  by  Slater, 


2  V  (d2pc/dV2\ 

3  2  (  dpc/dV  J 


(3) 


and  that  of  Dugdale  and  MacDonald, 

V  fd2pc/dV 2  -  10pc/9V2\ 
2  \  dpc/dV  +  2pc/3V  ) 


bear  a  notable  similarity  as  presented,  it  was  Vashchenko  and  Zubarev  (4)  who  were  able 
to  present  a  unified  view  of  these  results,  along  with  the  tree-volume  theory.  They 
present  a  single  parameterized  algebraic  expression  which  accurately  captures  the  forms 
of  all  three  of  these  models,  depending  upon  the  value  of  the  model's  parameter.  This 
expression  is  given  as 


T 


d 2 

4-3 n  V 


( PcVr 


6 


2 


(5) 
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where  n  is  the  model  parameter.  When  n  takes  on  values  of  0,  2/3,  and  4/3  respectively, 
the  models  of  Slater  (2),  Dugdale  and  MacDonald  (3),  and  the  free-volume  theory  (4)  are 
respectively  recovered. 

To  arrive  at  equation  5,  Vashchenko  and  Zubarev  postulated  a  hypothetical  form  for  the 
Poisson  ratio,  v,  in  terms  of  the  cold  pressure,  volume,  and  model  parameter  n: 

1  -  u0  (  <2  In  V  \ 

l  +  ^o  (  ndlnpc) 

1  +  ^oV  n  d\npc  ) 

Coupled  with  certain  relations  specifying  the  manner  in  which  the  Poisson  ratio  relates 
to  the  lattice  vibrational  frequencies,  equation  6  was  shown  to  lead  directly  to  equation  5. 

Therefore,  if  one  accepts  their  postulated  form  for  the  Poisson  ratio,  equation  6,  then  the 
three  Griineisen  models  of  interest  (2,  3,  4)  follow  from  equation  5  as  specific  instances  of 
the  proposed  form.  Pastine  ( 8 )  has  also  utilized  an  explanation  involving  the  influence  of 
the  Poisson  ratio  as  a  means  to  interpret  and  differentiate  expressions  for  the  Griineisen 
formulation. 

Conceptually,  such  an  approach  makes  sense,  and  there  is  certainly  nothing  invalid 
about  it.  However,  the  intuitive  understanding  that  one  may  obtain  from  equation  6  is, 
at  best,  quite  limited,  especially  when  the  Poisson  form  is,  itself,  expressed  in  terms  of 
volume  derivatives  of  the  cold  pressure. 

What  is  more,  the  original  derivations  of  the  models  of  Slater,  Dugdale  and  MacDonald, 
and  the  free-volume  theory  were  not  made  on  the  basis  of  this  or  any  other  Poisson  ratio 
formulation.  That  researchers  could,  after  the  fact,  discover  the  Poisson  ratio  as  a  means 
to  differentiate  the  models  was  a  useful  and  convenient  resolution  to  the  useful  question 
of  "what  scientific  explanation  can  explain  the  functional  difference  between  these  three 
models?"  But,  because  the  original  models  were  never  derived  in  terms  of  Poisson  ratio, 
the  Poisson  explanation  fails  to  answer  the  more  fundamental  question  of  "what  is  it 
about  the  formulation  of  these  three  models  that  leads  to  their  functional  differences?" 
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3.  Result 


In  this  report,  a  different  approach  is  taken  to  compare  these  historical  models,  based 
upon  treating  the  material  lattice  as  a  spring-like  structure,  in  which  the  various  models 
may  be  differentiated  on  the  basis  of  what  kind  of  spring  composes  the  lattice  and  what 
measure  is  taken  to  quantify  the  spring's  characteristic  spatial  dimension.  The  key  link 
in  eventually  relating  the  spring's  response  back  to  the  Griineisen  function  is  the 
classical  equation  2,  which  defines  the  Griineisen  function  in  terms  of  the  characteristic 
vibrational  frequency  of  the  lattice  u. 

Consider  a  lattice  with  spacing  A  such  that  A3  =  V .  When  atomic  vibration  is  not 
considered  (i.e.,  at  0°  temperature),  the  potential  energy  of  the  lattice  network  is  defined 
by  the  cold  energy  of  the  lattice  Ec,  which  will  be  a  function  of  only  a  spatial  dimension 
that  is  characteristic  of  the  lattice,  call  it  £.  This  characteristic  dimension  might  take  on 
units  of  length  (as  in  the  lattice  spacing  A),  though  it  need  not  necessarily  do  so;  rather,  it 
could  take  the  form  of  volume  V  or  some  other  spatially  related  quantity 

For  a  quantity  characterized  by  an  isotropic  field  potential,  such  as  the  lattice  energy  Ec, 
the  negative  of  the  gradient  of  the  potential  with  respect  to  £,  in  this  case  — Vi?c(£),  is  an 
indicative  measure  of  the  field's  "force,"  while  the  Laplacian  of  Ec,  that  is  V2i?c(£),  is  a 
measure  of  the  field's  "stiffness."  The  terms  "force"  and  "stiffness"  are  quoted  here 
because  depending  on  the  nature  of  the  characteristic  dimension  £,  the  physical  units  of 
the  characteristic  lattice  "force"  might  be  other  than  force  and  likewise  for  the  "stiffness." 

Vibration  analysis  indicates  that  the  system's  characteristic  stiffness  is  associated  with  a 
characteristic  vibrational  frequency  according  to  the  following  relation: 

V2£c(0  oc  cf  ,  (7) 

where  the  symbol  "oc"  denotes  a  direct  proportionality.  This  expression  has  been  used  in 
the  study  of  lattice  mechanics  (9, 10).  As  in  the  case  of  the  characteristic  force  and 
stiffness,  the  £-based  characteristic  frequency  c may  not  have  the  physical  units  of 
1/time  that  are  traditionally  associated  with  frequency,  though  it  will  have  dimensions  of 
propagation  speed  per  unit  characteristic  dimension.  However,  it  may  be  related  to  a 
characteristic  frequency  u  that  does  have  traditional  units  (i.e.,  1/time).  This  connection 
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is  made  through  the  definitional  relations  of  vibration  theory  which  connect  the 
characteristic  frequency  and  propagation  velocity  through  the  vibrational  wavelength: 


u£  =  C  =  u\  .  (8) 

Thus,  combining  equations  7  and  8  gives  the  following  proportion: 

V2£c(0  oc  (y)  .  (9) 

Equation  9  is  the  simple  relation  from  which  the  results  of  this  report  flow.  Applying  it 
to  the  task  of  acquiring  a  relationship  for  the  Griineisen  function  is  accomplished  by 
starting  with  equation  9  and  isolating  ca2.  Take  the  derivative  of  this  equation  with 
respect  to  volume  V  and  divide  the  resulting  equation  by  the  original  equation.  The 
result,  in  light  of  equation  2,  is  given  as 

r  v  W_  (vV^K) 

2 


There  are  really  two  parameters  that  characterize  equation  10.  The  first  is  the 
characteristic  lattice  dimension  £,  and  the  second  is  physical  nature  of  the  lattice  spring 
whose  stiffness  is  characterized  by  V2 /£.(£).  Depending  on  whether  the  lattice  spring  is 
considered  a  simple  one-dimensional  (1-D)  spring  or  a  three  dimensional  (3-D)  radial 
spring,  the  Laplacian  operator  would  be  expressed  as 

{d2Ec/d£2  (simple  1-D  spring) 

(11) 

d2Ec/d £2  +  (2 /£)dEc/d£  (3-D  radial  spring) 


From  equation  10,  in  light  of  equation  11,  it  may  be  directly  shown  that  Slater's  relation 
is  recovered  when  the  lattice  spring  is  a  simple  1-D  spring  and  the  characteristic  spatial 
dimension  £  of  the  lattice  is  taken  as  its  cell  volume  V.  Similarly,  the  relation  of  Dugdale 
and  MacDonald  is  recovered  by  taking  the  lattice  spring  to  be  simple  1-D,  while  taking 
the  characteristic  spatial  dimension  £  as  the  lattice  spacing  A.  Finally,  the  free-volume 


5 


theory  may  be  recovered*  from  equation  10  when  the  lattice  spring  is  taken  as  a  3-D 
radial  spring  and  the  characteristic  spatial  dimension  £  of  the  lattice  is  the  lattice  spacing 
A.  Table  1  categorizes  some  of  the  relevant  quantities  for  these  three  models. 


Table  1.  Characteristic  spring  type  and  spatial  dimension  for  various  Grtineisen  theories  and 
physical  units  of  associated  characteristic  quantities. 


Theory 

Lattice  Spring 
Type 

Characteristic 

Spatial  Dimension,  £ 

Units  of  Characteristic 

Force 

Stiffness 

Slater 

Simple  (1-D) 

V 

Pressure 

Pressure /volume 

Dugdale-MacDonald 

Simple  (1-D) 

A 

Force 

Force /length 

Free-volume  theory 

Radial  (3-D) 

A 

Force 

Force /length 

These  differentiating  characteristics  are  not  merely  coincidental,  which  just  happen  to 
produce  the  given  models;  rather,  they  are  indicative  of  the  logic  and  methodologies 
employed  by  the  researchers  in  their  original  derivations.  In  the  case  of  Slater's 
derivation,  a  constant  Poisson  ratio  was  assumed.  When  this  is  stipulated  for  a  material, 
the  elastic  and  bulk  moduli  remain  proportional.  This  allowed  the  vibrational  frequency 
to  be  cast  in  terms  of  volume  derivatives  of  pressure.  Thus,  it  is  not  surprising  that  the 
characteristic  dimension  of  the  Slater  model  is  volume  V,  leading  to  a  stiffness  expressed 
in  terms  of  pressure  per  volume. 

In  the  case  of  the  Dugdale-MacDonald  model,  the  authors'  explicit  contention  was  that 
the  characteristic  lattice  velocity  is  always  of  a  form  consistent  with  a  1-D  spring 
formulation  (which  is  to  say,  with  a  characteristic  stiffness  having  units  of  force  per 
length).  Plendl  (22)  has  also  asserted  this  relationship.  Dugdale  and  MacDonald  even 
went  so  far  as  to  show  that  Slater's  model  failed  to  give  the  proper  Grtineisen  value  of 
zero  for  a  harmonic  lattice  of  simple  1-D  springs.  Thus,  wholly  consistent  with  the 
model's  underlying  principal,  a  lattice  of  simply-connected  1-D  springs  would 
necessarily  retain  A  as  the  characteristic  dimension. 

The  free-volume  theory  differs  from  the  Dugdale-MacDonald  relation  only  in  that  the 
lattice  spring  is  taken  as  a  3-D  radial  spring  rather  than  as  a  simple  1-D  spring, 
reflecting  a  more  complex  view  of  reality  in  terms  of  the  interatomic  potential  as  a  3-D 
field  and  not  as  a  1-D  spring. 

*  Recognition  of  the  compatibility  of  the  free-volume  theory  with  the  underlying  tenets  of  equations  2 
and  7  was  also  noted  in  unpublished  work  by  Shanker  et  al.  (11). 
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It  is  important  to  note  that  when  £  defines  an  orientation-independent  lattice  measure 
such  as  A  or  V,  equation  7  describes  a  volumetric  oscillation,  since  an  ordinary  derivative 
of  Ec,  with  respect  to  the  characteristic  dimension  £,  implies  that  £  is  changed  between 
every  atom  in  the  lattice  network  (and  in  every  component  direction),  in  order  to 
evaluate  the  resultant  change  in  Ec.  Therefore,  in  such  cases,  the  quantity  may  be 
justifiably  termed  a  "volumetric  frequency,"  regardless  of  whether  the  springs  are 
characterized  as  1-D  or  3-D. 

However,  it  has  been  shown  (13)  that  when  the  lattice  interactions  are  restricted  to 
nearest  neighbors  only  (or  alternately  when  the  lattice  component  stiffness  is  purely 
harmonic  in  nature),  the  stiffness  associated  with  a  volumetric  change  is  in  proportion  to 
the  stiffness  associated  with  a  1-D  vibrational  wave  pulse.  Therefore,  under  the 
assumption  of  nearest-neighbor  interactions,  this  volumetric  frequency  is  equally 
indicative  of  a  traveling-wave  vibrational  frequency  from  which  the  T  function  actually 
arises. 

In  all  the  models  examined  in  this  report  (1,  3,  4),  the  original  authors  explicitly  invoked 
nearest-neighbor  arguments  in  their  derivations.  This  widely  used  nearest-neighbor 
assumption  (1,  3,  4, 14, 15)  serves  to  justify  the  use  of  in  equation  8  as  a  means  to  link 
the  vibrational  frequencies  of  volumetric  change  to  those  of  traveling  waves.  It  must  be 
remembered,  however,  that  in  the  general  case  when  distant  interactions  of  non-nearest 
neighbors  are  considered,  the  volumetric  frequency  and  vibrational  frequency  are  not 
identical  (13).  As  such,  equation  7  remains  merely  an  approximation  when  applied  to 
the  global  (as  opposed  to  the  pairwise)  potential. 


4.  Conclusion 


The  simple  direct  formulation  of  the  Griineisen  function  that  herein  arises  from 
elementary  vibration  analysis  serves  to  provide  a  clear  foundational  understanding  of 
three  popular  Griineisen  models  (1,  3,  4)  and  their  differentiating  characteristics.  Unlike 
the  historical  explanations  (4,  8)  that,  after  the  fact,  invoked  the  functional  behavior  of 
the  Poisson  ratio  as  the  means  to  functionally  differentiate  the  historical  models,  the 
current  approach  offers  an  explanation  in  terms  of  two  fundamental  modeling 
assumptions — the  nature  of  the  interatomic  "spring"  and  the  characteristic  spatial 
measure  of  the  lattice.  Such  a  differentiation,  by  its  nature,  provides  a  satisfactory 
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answer  to  the  question  of  "what  is  it  about  the  formulation  of  these  three  models  that 
leads  to  their  functional  differences?"  Such  understanding  is  absent  with  the  historical 
explanations,  where  the  differentiating  feature  is  cast  in  terms  of  the  Poisson  ratio. 
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D  LAMBERT 

EGLIN  AFB  FL  32542-6810 


20  SANDIA  NATL  LABS 

ATTN  MAIL  SERVICES  MS-0100 
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PROVIDENCE  RI  02912 

3  CALTECH 

J  SHEPHERD  MS  105  50 
A  P  INGERSOLL  MS  170  25 
T  J  AHRENS  MS  252  21 
1201  E  CALIFORNIA  BLVD 
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SCHOOL  OF  MATL  SCI  &  ENGNG 
K  LOGAN 
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KURZZEITDYNAMIK 
ERNST  MACH  INSTITUT 
H  ROTHENHAEUSLER 
H  SENF 

E  STRASSBURGER 
KLINGELBERG  1 
D79588  EFRINGEN-KIRCHEN 
GERMANY 


3  FRENCH  GERMAN  RSRCH  INST 
G  WEIHRAUCH 
R  HUNKLER 
E  WOLLMANN 
POSTFACH  1260 
WEIL  AM  RHEIN  D-79574 
GERMANY 

2  IABG 

M  BORRMANN 

H  G  DORSCH 

EIN STEIN STRA S SE  20 

D  8012  OTTOBRUN  B  MUENCHEN 

GERMANY 

1  INGENIEURBUERO  DEISENROTH 
AUF  DE  HARDT  33  35 
D5204  LOHMAR  1 
GERMANY 

1  NORDMETALL  GMBH 
L  W  MEYER 
EIBENBERG 
EINSIEDLER  STR  18  H 
D-09235  BURKHARDTSDORF 
GERMANY 

2  TU  CHEMNITZ 

L  W  MEYER  (X2) 

FAKULTAET  FUER  MASCH1NENBAU 

LEHRSTUHL  WERKSTOFFE  DES 

MASCHINENBAUS 

D-09 107  CHEMNITZ 

GERMANY 

1  TU  MUENCHEN 
E  IGENBERGS 
ARCISSTRASSE  21 
8000  MUENCHEN  2 
GERMANY 

2  UNIVERSITAET  PADERBORN 
DEPARTMENT  PHYSIK 

O  SCHULTE 
W  B  HOLZAPFEL 
D-33095  PADERBORN 
GERMANY 

1  BHABHA  ATOMIC  RSRCH 
CENTRE 

HIGH  PRESSURE  PHYSICS  DIV 
N  SURESH 

TROMBAY  BOMBAY  400  085 
INDIA 


19 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


1  NATIONAL  GEOPHYSICAL 
RSRCH  INST 
G  PARTHASARATHY 
HYDERABAD-500  007  (A.  P.) 

INDIA 

1  UNIV  OF  ROORKEE 

DEPT  OF  PHYSICS 
NDASS 

ROORKEE-247  667 
INDIA 

5  RAFAEL  BALLISTICS  CTR 
E  DEKEL 

Y  PARTOM 

G  ROSENBERG 
Z  ROSENBERG 

Y  YESHURUN 
PO  BOX  2250 
HAIFA  31021 
ISRAEL 

1  SOREQ  NUCLEAR  RSRCH 
CENTRE 

ISRAEL  ATOMIC  ENERGY 

COMMISSION 

Z  JAEGER 

81800  YAVNE 

ISRAEL 

1  ESTEC  CS 

D  CASWELL 
BOX  200  NOORDWIJK 
2200  AG 
NETHERLANDS 

2  EUROPEAN  SPACE  AGENCY  ESTEC 
L  BERTHOUD 

M  LAMBERT 

POSTBUS  BOX  299  NOORDWIJK 

NL2200  AG 

NETHERLANDS 

4  PRINS  MAURITS  LAB 
H  J  REITSMA 
E  VAN  RIET 
H  PASMAN 
R  YSSELSTEIN 
TNO  BOX  45 
RIJSWIJK  2280AA 
NETHERLANDS 


1  ROYAL  NETHERLANDS  ARMY 
J  HOENEVELD 

V  D  BURCHLAAN  3 1 
PO  BOX  90822 

2509  LS  THE  HAGUE 
NETHERLANDS 

1  INST  OF  PHYSICS 

SILESIAN  TECHNICAL  UNIV 

E  SOCZKIEWICZ 

44-100  GLIWICE 

UL.  KRZYWOUSTEGO  2 

POLAND 

1  INST  OF  CHEMICAL  PHYSICS 

A  YU  DOLGOBORODOV 
KOSYGIN  ST  4  V  334 
MOSCOW 
RUSSIA 

4  INST  OF  CHEMICAL  PHYSICS 

RUSSIAN  ACADEMY  OF  SCIENCES 
G  I  KANEL 
A  M  MOLODETS 
S  V  RAZORENOV 
A  V  UTKIN 

142432  CHERNOGOLOVKA 
MOSCOW  REGION 
RUSSIA 

1  INST  OF  EARTHS  CRUST 
P  I DOROGOKUPETS 
664033  IRKUTSK 
RUSSIA 

3  INST  OF  MECH  ENGINEERING 
PROBLEMS 

V  BULATOV 
D  INDEITSEV 

V  MESCHERYAKOV 
BOLSHOY,  61,  V.O. 

ST  PETERSBURG  199178 
RUSSIA 

1  INST  OF  MINEROLOGY  & 
PETROGRAPHY 

V  A  DREBUSHCHAK 
UNIVERSITETSKI  PROSPEKT,  3 
630090  NOVOSIBIRSK 
RUSSIA 
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2  IOFFE  PHYSICO  TECHNICAL 
INST  DENSE  PLASMA 
DYNAMICS  LAB 
E  M  DROBYSHEVSKI 
A  KOZHUSHKO 
ST  PETERSBURG  194021 
RUSSIA 

1  IPE  RAS 

A  A  BOGOMAZ 
DVORTSOV  AIA  NAB  18 
ST  PETERSBURG 
RUSSIA 

2  LAVRENTYEV  INST 
HYDRODYNAMICS 
L  A  MERZHIEVSKY 
V  V  SILVESTROV 
630090  NOVOSIBIRSK 
RUSSIA 

1  MOSCOW  INST  OF  PHYSICS  &  TECH 
S  V  UT YU  ZHNIKO V 
DEPT  OF  COMPTNL  MATHEMATICS 
DOLGOPRUDNY  1471700 
RUSSIA 

1  RESEARCH  INST  OF  MECHANICS 
NIZHNIY  NOVGOROD  STATE  UNIV 
A  SADYRIN 

P  R  GAY  ARINA  23  KORP  6 
NIZHNIY  NOVGOROD  603600 
RUSSIA 

2  RUSSIAN  FEDERAL  NUCLEAR 
CTR  -  VNIIEF 

L  F  GUDARENKO 
R  F  TRUNIN 
MIRA  AVE  37 
SAROV  607190 
RUSSIA 

1  ST  PETERSBURG  STATE  TECHNICAL 
UNIV 

FACULTY  OF  PHYS  AND  MECHS 
DEPT  OF  THEORETICAL  MECHS 
ATTN  AM  KRIVTSOV 
POLITECHNICHESKAY A  STREET  29 
195251  ST-PETERSBURG 
RUSSIA 


1  SAMARA  STATE  AEROSPACE  UNIV 
L  G  LUKASHEV 
SAMARA 
RUSSIA 

1  TOMSK  STATE  UNIV 
A  G  GERASIMOV 
5-TH  ARMY  ST  29-61 
TOMSK  634024 
RUSSIA 

5  DEPARTMENTO  DE  QUIMICA 

FISICA  I  FACULTAD  DE  CIENCIAS 
QUIMICAS 

UNIVERSIDAD  COMPLUTENSE  DE 

MADRID 

V  G  BAONZA 

M  TARAVILLO 

J  E  F  RUBIO 

J  NUNEZ 

M  CACERES 

28040  MADRID 

SPAIN 

1  UNIVERSIDAD  DE  CANTABRIA 
FACULTAD  DE  CIENCIAS 
DEPARTMENTO  DE  FISICA 
APLICADA 
J  AMOROS 

AVDA  DE  LOS  CASTROS  S/N 
39005  SANTANDER 
SPAIN 

1  CARLOS  III  UNIV  OF  MADRID 
C  NAVARRO 
ESCUELA  POLITEENICA 
SUPERIOR 
C/.  BUTARQUE  15 
2891 1  LEGANES  MADRID 
SPAIN 

1  UNIVERSIDAD  DE  OVIEDO 
FACULTAD  DE  QUIMICA 
DEPARTMENTO  DE  QUIMICA 
FISICA  Y  ANALITICA 
E  FRANCISCO 

AVENIDA  JULIAN  CLAVERIA  S/N 

33006  -  OVIEDO 

SPAIN 
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1  DYNAMEC  RSRCH  AB 
A  PERSSON 

PO  BOX  201 

S-15123  SODERTALJE 

SWEDEN 

7  FOI 

SWEDISH  DEFENCE  RSRCH  AGENCY 

GRINDSJON  RSRCH  CENTRE 

L  GUNNAR  OLSSON 

B  JANZON 

G  WIJK 

R  HOLMLIN 

C  LAMNEVIK 

L  FAST 

M  JACOB 

SE- 14725  TUMBA 

SWEDEN 

2  SWEDISH  DEFENCE  RSRCH  ESTAB 
DIV  OF  MATERIALS 

S  J  SAVAGE 
J  ERIKSON 
S-17290  STOCKHOLM 
SWEDEN 

2  AWE 

M  GERMAN 
W  HARRISON 
FOULNESS  ESSEX  SS3  9XE 
UNITED  KINGDOM 

5  DERA 
I  CULL1S 
J  P  CURTIS  Q13 
AHARTQ13 
K  COWAN  Q13 
M  FIRTH  R31 
FORT  HALSTEAD 
SEVENOAKS  KENT  TNI 4  7BP 
UNITED  KINGDOM 


1  K  TSEMBELIS 

SHOCK  PHYSICS  GROUP 
CAVENDISH  LAB 

PHYSICS  &  CHEMISTRY  OF  SOLIDS 
UNIV  OF  CAMBRIDGE 
CAMBRIDGE  CB3  OHE 
UNITED  KINGDOM 

1  L  VOCADLO 

DEPT  EARTH  SCIENCES 
UNIV  COLLEGE  LONDON 
GOWER  ST 
LONDON  WC1E  6BT 
UNITED  KINGDOM 

7  INST  FOR  PROBLEMS  IN 

MATERIALS  SCIENCE 
S  FIRSTOV 
B  GALANOV 
O  GRIGORIEV 

V  KARTUZOV 

V  KOVTUN 

V  MILMAN 

V  TREFILOV 

3,  KRHYZHANOVSKY  STR 
252142  KIEV- 142 
UKRAINE 

1  INST  FOR  PROBLEMS  OF  STRENGTH 
G  STEPANOV 

TIMIRY  AZE  V  SKA  YU  STR  2 

252014  KIEV 

UKRAINE 


1  UK  MINISTRY  OF  DEFENCE 
G  J  CAMBRAY 

CBDE  PORTON  DOWN  SALISBURY 
WITTSHIRE  SPR  OJQ 
UNITED  KINGDOM 
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